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Abstract. - We study the steady-state response to applied stress in a simple scalar model of 
sheared colloids. Our model is based on a schematic (F2) model of the glass transition, with 
a memory term that depends on both stress and shear rate. For suitable parameters, we find 
transitions from a fluid to a nonergodic, jammed state, showing zero flow rate in an interval of 
applied stress. Although the jammed state is a glass, we predict that jamming transitions have 
an analytical structure distinct from that of the conventional mode coupling glass transition. 
The static jamming transition we discuss is also distinct from hydrodynamic shear thickening. 



When subjected to an applied shear stress, concentrated hard-sphere colloids are known to 
shear-thicken [20010] . An understanding of this in terms of hydrodynamic interactions has 
been developed jolGj . Recently, a different (although seemingly related) phenomenon has been 
found. Bertrand and co-workers [7] report an experiment in which shearing a concentrated 
suspension induces a transition from a fluid to a metastable solid which persists after cessation 
of shear. Vibration of the sample restores it to a fluid state. 

Because the solid paste persists in the absence of flow, such 'static jamming' cannot be 
explained by hydrodynamic forces alone. An alternative view is that after cessation of flow 
the material remains stressed; this stress is what maintains the arrest. Indeed, jammed 
suspensions often take on a lumpy dry appearance as particles protrude partially from the 
fluid surface (a dilatancy effect). This entails large capillary stresses which might maintain 
the jam; lumps of paste removed from a rheometer would then remain solid. Vibration 
could destroy the jammed structure and restore both the fluidity and a wet appearance, as 
is observed [3121121 • Rather than model all this directly, we focus here on the simpler case of 
stress-induced jamming within an idealised rheometer, with no free surfaces. 

It was recently proposed that jamming in granular materials ^U] might be closely connected 
to the glass transition |11U12| . If so, the jammed colloidal state found by Bertrand et al could 
resemble a glass, albeit one which is stressed and hence anisotropic. (Colloidal glasses are 
well-documented A speculative scenario for the formation of a shear-induced glass is as 

follows: the applied shear stress alters the structure of the material, initially through shear flow 
and hydrodynamics. But as the material thickens, its structure becomes more susceptible to 
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dynamical arrest, because an increased number of close interparticle contacts hinders diffusion. 
This promotes a nonergodicity transition, at concentrations somewhat lower than would be 
required without the stress. In the arrested state, flow ceases but stress remains; the stress is 
now sustained by interparticle and/or entropic forces rather than hydrodynamics. 

The consideration of such a phenomenon is complicated by the tensorial nature of applied 
stresses; a system jammed by one component of an applied stress might be refluidised by 
another |1U| . In the current work, we ignore these tensorial complications, and treat jamming 
as a generalised glass transition within a scalar model. (For a somewhat complementary 
approach, see This model builds on the ideas of the mode-coupling theory (MCT) of 

the glass transition, but introduces two new features: stress-induced arrest, and strain-induced 
memory loss. The combination of these features allows new shear-thickening and jamming 
scenarios to emerge. Fuller details of our calculations will appear elsewhere |15j . 

For a comprehensive description of MCT, see We outline only those aspects pertinent 
to the present work. The central quantities of the theory are density fluctuations at wavevector 
q, Sp(q,t). The correlators of these quantities, q (i) = (5p(q, t)Sp(— q, 0))/(|6p(q)| 2 ) may be 
measured in scattering experiments, and provide a description of the system's dynamics. In 
a liquid, the system is ergodic and </> q (i) decays to zero with time for all q. In a glass it 
does not: lim^oo 4>q(t) — / q > 0, where / q is a nonergodicity parameter characteristic of the 
arrested structure. A finite / q represents the inability of the structure to relax on lengthscale 
~ 27r/g, preventing an initial fluctuation from fully decaying. 

Equations of motion can be found for the correlators q (i). On making approximations 
suitable for a colloidal system, and dropping q-subscripts, the result is [T7|: 



where r Q sets the timescale for the microscopic dynamics. Here m(t — t') is the memory func- 
tion, and describes a retarded friction effect which, in the colloidal glass transition, arises by 
caging of a particle by its neighbours. In MCT, the memory function is found approximately 
by integrating (over wavevectors) a quadratic product of correlators, with coupling constants 
that depend on the static structure factor S(q) — (1/N){\5p(q)\ 2 ) of the system. This ap- 
proximation means that anharmonic interactions between density fluctuations are ignored. 

Generalisation of the wavevector-dependent MCT to sheared systems is challenging, and 
so far has shown only shear-thinning behavior 18, 19,20 . Analogous behaviour is seen in 
mean field spin models of driven glasses [211 . To address the more complicated issue of jam- 
ming, a simpler starting point is needed. A promising one lies with 'schematic models' which 
capture many key features of the MCT glass transition. These schematic models consider a 
single correlator, rather than the infinite set {4>q{t)}. The memory function is then written 
as a polynomial of the correlator, with coefficients (coupling constants) that schematically 
represent system variables. Larger coupling constants correspond to higher densities or lower 
temperatures. We choose to build upon a simple schematic model (the 'F2 model') for which 
m(t) — v<j) 2 (t). This model has an idealised glass transition at v = 4 [22]: for v < 4, the only 
admissable solutions have / = 0, while for larger v nonergodic solutions, with / > 0, appear. 

To incorporate shearing into this picture we need to account for two effects. Firstly, flow 
suppresses memory: particles separated by a small distance in the flow gradient direction be- 
come well separated for times beyond I/7 (with 7 the shear rate). Beyond this timescale, any 
system at nonzero flow rate should lose memory of previous configurations, so that ergodicity 
is restored |23II21| . (In q-space, the important fluctuations, which have wavevectors near the 
peak of S(q), are advected on this timescale to higher q where they decay rapidly. Hence cages 
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are destroyed, and even those fluctuations that are not directly advected, become ergodic |18|.) 
The second effect is jamming: the fact that stress (as distinct from flow) can promote arrest. 
Stress can hinder diffusion by distorting cages and creating more close contacts between par- 
ticles, particularly if these tend to arrange into a load-bearing structure JO]- (This idea does 
not exclude a strong transient role for hydrodynamic effects in determining what structures 
actually do arise.) Such physics is unlikely to be captured within the harmonic approximation 
normally used in full MCT. To capture it, one might make a harmonic expansion using the 
S(q) of some anisotropic reference state in which the thermodynamic stress was nonzero: at 
the schematic level, we should thus allow a stress dependence of the coefficients, alongside an 
explicit memory loss arising from nonzero 7. 

To allow for both these effects, we model the memory function under shear as follows: 



Here the parameter v represents the system's tendency to arrest in the absence of any external 
forcing, and a represents the degree to which this intrinsic memory is enhanced by a shear 
stress of magnitude a. (Clearly the interesting case for jamming systems is a > 0.) At small er 
one could argue for a quadratic dependence (since positive and negative shear stresses should 
be equivalent) but we shall find that the interesting behavior of the model in fact lies in a 
window of stress where a is not small. (The chosen a dependence can then be interpreted 
as a linearisation about some point near the middle of this window.) The exponential form 
chosen for flow-induced memory loss is somewhat arbitary, but the results do not depend 
much upon this choice: the important part of this function involves jt > 1 (so that the 
apparent nonanalytic dependence at small flow rates is again immaterial). Another form, 
m(t) oc 1/(1 + {jt) 2 ), which is better motivated by microscopic considerations yields 
qualitatively similar numerical results, but the exponential is preferable in analytic work |15|. 

In steady flow the stress a and strain rate 7 are not independent, but are related through 
a viscosity 77(7) = er/7 which depends on the state of the material. To close our schematic 
model, we need a prescription for this viscosity. In linear response, the viscosity may be 
expressed as the time integral of a stress correlator |24j . In our simplified model, we have only 
one correlator, and so write 



where we choose units so that r\ equates to a characteristic relaxation time r. Equation is 
again somewhat ad-hoc, since we are not studying the linear response regime; but equating 
77 to an integral of 4>(t) 2 (say) would give very similar results. More generally, according to 
the full MCT the viscosity diverges like the relaxation time of a typical correlator as one 
approaches a colloidal glass transition (a correspondence that is not lost under shear |18j): 
Eq[21 captures this. It means that dynamical arrest (nonzero /) implies a divergent viscosity 
and hence zero shear rate. This is fully consistent with arguments made above that finite 7 
prohibits arrest. (Note that the shear rate refers to the steady state value, so this does not 
exclude the possibility of sublinear creep of the strain, with vanishing shear rate at late times.) 

Our model is completely defined by Eqs. El and It has been solved numerically by 
adapting an established algorithm for MCT equations |2SJ. This allows iteration (from above 
or below) of the relaxation time r (or equivalently 7) to a self-consistent solution for fixed 
v a , a and a. Below, when we refer to stability of the solutions, we mean stability under this 
iteration. (Physical stability is discussed at the end.) The numerical results were also checked 
by analysis of limiting cases |15|. Note that unstable solutions (in the sense just defined) also 
exist: by inspection, for 7 = 0, the model reduces to the standard F2 model, with coupling 



m{t) = (v + aa) exp(— jt)^ 2 ^). 
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Fig. 1 - Flow curves, stress a versus shear rate 7, for a = 0.95. For the two largest values of v , it 
appears that for a window in a, the relaxation time has diverged. Analytic calculations of the limits 
of this window are indicated as horizontal lines near the stress axis. These values of the stress are 
dubbed o c \ and <r C 2, as shown here for one of the parameter sets. Larger shear rates (which have been 
made dimensionless by mutiplication with the 'bare' relaxation time r D ) are found at larger stresses. 



v — v a + aa. Therefore, for nonzero a, a nonergodic solution at zero shear rate always formally 
exists for a sufficiently large (v + ever > 4) . Numerically however, it was found that this state 
is unstable with respect to an ergodic solution whenever the latter exists. However, for some 
parameter choices, there is no ergodic solution within a certain window of stress. Here the 
nonergodic jammed state, with finite stress but zero flow rate, is stable. 

The resulting 'full jamming' scenario (see also |14) 1 forms part of a wider range of rheolog- 
ical behavior. Fig. ^ shows three thickening scenarios, dependent on model parameters; v is 
chosen close to the quiescent glass transition and a is chosen close to unity. Upon increasing 
v a there is a progression from a monotonic, continuously shear-thickening curve, via a non- 
monotonic S-shaped curve, to a curve that extends right back to the vertical axis. It is known 
that any flow curve with negative slope is unstable to shear banding |2f>| ; hence for the latter 
two scenarios, in any experiment at controlled shear rate, the stress would be expected to 
jump discontinuously from the lower to the upper branch before reaching the point of infinite 
slope. Such discontinuous shear thickening is widely reported |1I2I3I4| and usually attributed 
to hydrodynamic interactions Our work suggests that, at least in some systems, this may 
not be the only mechanism at work. In particular, Fig. ^ shows shear thickening at Peclet 
numbers 7T0 ~ 10~ 4 , rather than values of order ten predicted by hydrodynamics. We will 
discuss this in greater depth in a future paper |15j . 

For the largest values of the parameter v Q , in Fig. ^ there is a range of stress for which the 
shear rate returns to zero: there is then no ergodic solution, and the jammed state is stable. 
This represents full, static jamming. The lower and upper endpoints cr c i and a c i of the stable 
jammed state represent distinct jamming transitions. Their critical stresses obey 



f c [(v + aa c ) f c - 2] = er c , 



(4) 
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Fig. 2 - 'Phase diagrams' for the model for various v . The lines denote transitions in the (a, cr)-plane. 
All states below the curve for a given value of v are fluid states, whilst those above (and on) the line 
are nonergodic, jammed states. 



where f c is given by the largest solution of jzrf — i v o + ota c )fc ■ Such transitions exist 
provided that both v D and a are sufficiently large. The meaning of v a is clear: it represents the 
tendency of an unsheared suspension to arrest, and is controlled by the colloid concentration 
and interactions jlfij . The experimental meaning of a is less clear. Bertrand et al [Jj found 
that, for concentrations below a certain value, their samples showed ordinary thickening, 
whilst above this value the shear-induced solid was seen. The behaviour illustrated in Fig. 
is reminiscent of this. Note that the re-fluidisation under increasing stress is dependent upon 
the value of a: if this parameter is sufficiently large (for a given value of v a ) this re-fluidisation 
is not present. This is illustrated in a 'phase diagram' of the model, dividing parameter space 
into ergodic and nonergodic regions (Fig. [2J. This shows that at large enough stresses, jammed 
(fluid) states arise for a > 1 (a < 1). However, for particle densities close to but below the 
quiescent glass transition, for a < 1 the system jams in an intermediate window of stress. 

We now examine the generic behavior at the jamming transitions within our model, and 
compare this to the glass transition as conventionally addressed in MCT ^^J. Note that, 
because of an expected shear-banding instability associated with the decreasing part of the 
flow curve |26| . the lower jamming transition at cr c \ might be more complicated to access 
experimentally than the upper one at a C 2- The latter should always be accessible by raising 
the stress in the fluid until discontinuous shear thickening occurs (taking the system to the 
upper branch of the flow curve) and then reducing the stress through the transition at a C 2- 
This practical difference aside, we find similar behavior at both transitions for all the features 
discussed below. 

The standard MCT glass transition is associated with a bifurcation of / in parameter 
space: at sub-critical coupling, there is no solution with a nonzero /. At the transition point, 
there is a bifurcation in the solution adopted by the dynamics for t — > oo. This leads to a 
jump in /, from zero to a finite value, followed by a nonanalytic variation as the coupling is 
increased. At our jamming transitions, there is no bifurcation in / at the transition point, and 
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Fig. 3 - Schematic variation of 7 and / on varying a. Solid (dashed) lines denote stable (unstable) 
solutions. Arrows indicate flow under iteration. Note that the ergodic, flowing solution remains stable 
beyond the bifurcation point where a nonergodic solution first arises: le, v + aa c i > 4. Hence there 
are no nonanalyticities after the discontinuity in / at the lower or upper jamming transition. 

so the jump in / is followed by a linear variation with stress. Indeed, the jamming transitions 
in our model mark the disappearance of the ergodic (flowing) solution - stable whenever 
it exists - rather than the first appearance of the nonergodic (arrested) one. The latter is 
where the bifurcation occurs, and while in the MCT of conventional glasses this coincides 
with disappearance of the ergodic solution, here it does not. The bifurcation and stability 
behavior of our model is summarised in Figure^ There is a second, closely related, difference 
between the jamming transitions described here and the conventional MCT glass transition. 
In the vicinity of a conventional transition, there are two divergent timescales, those of the 
(short time) (3 and (long time) a relaxation processes. In contrast, at the jamming transitions 
described here, there is only one. Again, this is because the transition occurs at a different 
point in parameter space to the underlying glass transition. 

Throughout the above discussion, we have tacitly assumed that the stability of our solu- 
tions under iteration also governs the physical selection mechanism: a flowing state is preferred 
to a jammed one, whenever both exist. Physically this is not obvious, but plausible; although 
the iteration does not map directly onto a dynamical evolution, it does suggest that any tran- 
sient violation of Eq0 leading temporarily to an infinitesimal shear rate in the nonergodic 
state, will carry the system towards a steadily flowing solution unless a c \ < a < a C 2- Nonethe- 
less, recall that a formal nonergodic solution does exist at zero shear rate whenever v +a<r > 4. 
Were this to be physically stable and the flowing solution unstable, we would have only one 
jamming transition, showing all the behaviour of a conventional glass transition. But reports 
of flowing colloids at high stresses in discontinuously shear-thickening systems make 
that outcome unlikely. A more plausible (though equally speculative) scenario is hysteresis: 
it might be posssible to maintain a metastable jammed state outside the range o c \ < a < cr c2 
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by starting within that range and carefully changing the stress. The model predicts no limit 
to this at high stresses, but at low ones, metastability cannot extend below a C Q = (4 — v a )/a. 

In this Letter we have presented a simple model of jamming in suspensions, based on a 
schematic mode-coupling model of the glass transition. Our model exhibits transitions to a 
nonergodic state, which has properties similar to that of the idealised glass state of MCT |lfij . 
but (assuming that the iterative stability of solutions also controls the physical behavior) the 
jamming transitions that we find have a novel structure. They are not controlled directly by 
a bifurcation of the nonergodicity parameter, but instead by the disappearance of an ergodic, 
fluid phase within a window of stress. Aspects of the model appear consistent with recent 
experiments 0, although more work is required to clarify the nature of static jamming in 
dense colloids under stress. We hope that our work will stimulate such experiments. 
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